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Abstract 

We show that domain walls are probes that enable one to distinguish large- 
distance modified gravity from general relativity (GR) at short distances. For 
example, low-tension domain walls are stealth in modified gravity, while they 
do produce global gravitational effects in GR. We demonstrate this by finding 
exact solutions for various domain walls in the DGP model. A wall with 
tension lower than the fundamental Planck scale does not inflate and has 
no gravitational effects on a 4D observer, since its 4D tension is completely 
screened by gravity itself. We argue that this feature remains valid in a generic 
class of models of infrared modified gravity. As a byproduct, we obtain exact 
solutions for super-massive codimension-2 branes. 
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1 Introduction and Summary 



A model of large distance modification of gravity pQ is based on the existence of 
a new mega-scale in the theory, r c ~ Hq 1 ~ 10 28 cm. In this model, the gravi- 
ton propagates extra degrees of freedom (two helicity-1 and one helicity-0 states), 
that lead to the van Dam-Veltman-Zakharov (vDVZ) [2] type discontinuity in the 
linearized approximation. 

However, the continuity in recovering the Einsteinian metric for localized sources 
in the r c — > oo limit, is due to non-linear interactions of these extra states, that 
become strongly coupled, and 'self-shield' at short distances [3]. 

Because of these non-linearities, any spherically-symmetric source of a Schwarzs- 
child radius r g ^ r c , gets endowed with a new physical distance scale r* ~ (r^r 2 ) 1 / 3 ^ 
r c . The modifications to the metric now depend on whether the source is bigger or 
smaller than tv For sources localized well within their own r*, the metric inside 
the ?vsphere (r <C r*) is almost Einsteinian, with tiny corrections, which can be 
computed either in the l/r c -expansion [HEBE] or exactly [7J. 

For the existing realistic localized sources, such as stars, or planets (or even 
galaxies) the values of r* are huge, and correspondingly, the deviations from the 
Einsteinian metric are strongly suppressed near the sources. One expects that the 
corrections are potentially measurable [5], [U [9j [7J by Lunar laser ranging experiments 
jlOj , but in general, detection of the modifications at short scales is extremely hard. 
In this respect, it is important to identify other types of gravitating sources, for 
which modifications are significant, even at the short distances. 

We shall show that domain walls, are exactly these type of sources. In particular, 
we will find that low tension walls do not gravitate in the DGP model pQ. We will 
also argue that the latter property should persist in other theories of large distance 
modified gravity. 

Another important question concerns the form of the metric beyond r*. The 
naive perturbative expansion suggests that for r* <C r <C r c the metric should 
have an approximately four-dimensional, 1/r, scalar-tensor-gravity type form [HE]. 
However, a non-perturbative solution of Ref. [7J exhibits a different behavior: for 
t ^> r* the metric turns into the one produced by a five dimensional source. This is 
because the four-dimensional mass of the source is completely screened by a halo of 
non-zero curvature that surrounds the source and extends to the distances ~ r* [7J. 

It is important to understand the significance and universality of the screening 
mechanism [7J. The question whether the direct crossover, beyond r*, into the five 
dimensional regime, is a general property of well-localized sources, should be under- 
stood better. Although, in the present work we won't be able to clarify this issue 
completely, interestingly enough, we will find a similar screening mechanism for the 
domain walls. As a result, properties of the domain walls are dramatically differ- 
ent in modified gravity, and they could serve as very effective probes to distinguish 
between modified gravity and general relativity (GR). 

Consider 4D space-time labelled by coordinates (t, x, y, z) and introduce a Nambu- 
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Goto domain wall with the stress-tensor 

= a 5(z) diag(l, —1, —1, 0) , (1.1) 

where a is the tension of the wall and z denotes the coordinate transverse to its 
worldvolume. For simplicity, we consider the wall to be straight in this coordinate 
system and localized at the point z = 0. For the moment we ignore the transverse 
size of the wall. Gravitational field of such a source was found by Vilenkin [llj and 
by Ipser and Sikivie (12] (VIS), and reads as follows: 

ds 2 = (1 - H\z\) 2 {-dt 2 + e 2Ht (dx 2 + dy 2 )) + dz 2 , (1.2) 

where H = 2i:G^a. One simple consequence of this metric is that a freely falling 
observer will be repelled from the wall. In general, an observer outside the wall floats 
in the Rindler space while on the wall he/she observes a 3D de Sitter expansion. 

Using the above knowledge, and trying to guess the behaviour of the domain 
walls in DGP, one comes to the following puzzle. Naive intuition would tell us 
that, at sufficiently short distances, the wall on the brane should behave as in 4D, 
and thus, inflate. On the other hand, at large distances, the wall should behave as 
a codimension-two object, which are known to produce a static metric with just a 
deficit angle, at least for sufficiently small tensions. How could a non-singular metric 
interpolate between the inflating and static patches? 

As we shall see, the resolution of this puzzle lies in the dramatically different 
behaviour of walls in GR and in DGP-like theories. In modified gravity such a 
domain wall becomes stealth, as long as its tension does not exceed a certain critical 
value. The wall will have no gravitational effects whatsoever. The presence of such 
a wall can only be detected by scattering on its core. For instance, if gravity is 
modified at cosmological distances, there could exist a domain wall passing through 
the Solar system; we would be able to discover it only by coming into a contact with 
its core. On the other hand, a long-range field of such a domain wall would be felt 
in the Solar system if gravity is described by GR. 

Detailed explanations of why this takes place is given in the bulk of the paper. 
Here we discuss one way to understand this in the context of DGP. Consider the 
modified Einstein equation: 

G<5 + m c {K^ - 9lxu K) = BkGnT^ , (1.3) 

where m c = r^T 1 , Gjf} is a 4D Einstein tensor, is an extrinsic curvature tensor 
of 4D space-time, and K is its trace (for the full set of equations see, e.g., [7]). The 
second term on the l.h.s. of Eq. (11.31) is what distinguishes it from the conventional 
Einstein equation. Moreover, (jl.3p can be interpreted as a GR equation with an 
effective source = T Mi , — 8 ^ (K^ u — g^K). For the domain wall the extrinsic 
curvature terms exactly compensate the stress tensor (II -ip giving rise to = 0. 
Hence, the tension of the wall, as seen from the point of view of a 4D observer, is 
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screened entirely by gravitational effects encoded in the extrinsic curvature. Not 
surprisingly, the domain wall worldvolume remains flat, and so does the metric on 
the brane: 

ds 2 \ 4 D = -dt 2 + dx 2 + dy 2 + dx 2 . (1.4) 

Furthermore, we shall see that this screening takes place inside the core of the 
wall. Hence, the analog notion of the r* scale for a domain wall (understood as 
where the self- shielding takes place) coincides with its thickness, 

r (dw) _ j 

1 * "core • 

This is to be compared to the Schwarzchild-like case [7J, where the shielding also 
occurs, and r* extends outside the source. The net result is the screening of the 4D 
tension/mass in both cases. 

Although the above arguments were given for the DGP model, as we show in 
Section HI a similar mechanism is expected to remain valid generically in models of 
long distance modified gravity in which the graviton is no longer massless [13j and 
acquires extra polarizations. 

2 Domain Walls in Perturbation Theory 

In this subsection we analyze the domain wall metric in perturbation theory. Sub- 
tleties of such an approach do manifest themselves already in the case of GR, as has 
been discussed in [T3]. Here, we reiterate and generalize those arguments. Consider 
again 4D space-time and a domain wall in it. Its metric is given by (11. 2p . Can this 
metric be recovered in a perturbative calculation? This question was addressed in 
Ref. [2] and we shall briefly review the argument for clarity. 

Given that the equation for the linearized metric is second order in derivatives, 
and from the symmetries of the source (II. ip . one naively obtains the result that 
the linearized metric potentials behave like oc \z\. In appropriate coordinates, the 
metric takes the form 

ds 2 \ Pert . Th. ^ (1 - 2H\z\) (-dt 2 + dx 2 + dy 2 ) + dz 2 , (2.1) 

where H = a/4m 2 p . This is consistent with the linearizatioin of the exact solution 

ds 2 \ U n. ~ (1 - 2H\z\) (-dt 2 + (1 + 2Ht)(dx 2 + dy 2 )) + dz 2 . (2.2) 

Indeed, one can check that (12.11) and (12.21) differ by a gauge transformation d^ u ) 
with £^ of the form £ 2 = 0, £° = H(x 2 + y 2 )/2 and = Htx % (here x % label x 
and Hence, to linear order in perturbation theory, the metric does not capture 
the fact that the wall inflates, which makes sense since the curvature scalar of the 

5 Note, though, that this transformation diverges at infinity. Hence, (conserved) sources that 
do not decay at infinity can probe the difference between (|2.ip and (|2.2[) . 
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worldvolume in the full solution is of order H 2 . The general lesson to be learned 
from the above considerations is as follows: it seems that the linearized theory does 
not capture all the physics, and one needs to go to higher orders in perturbation 
theory (or possibly to the full nonlinear theory). 

Now, let us perform the perturbative calculations for a domain wall localized 
on the brane in the DGP model. Here, we consider the walls with a sub-critical 
tension a <^ M 3 . The subtleties outlined above manifest themselves in this case 
too. The linearized equations have a family of solutions, some of them are static 
and some are time dependent. It is hard to guess a priori what the right linearized 
solution is. However, we know from the exact solution of the next section the 
following: Unlike the domain wall in 4D GR, the worldvolume of a domain wall 
with cr ^ M'l in DGP does not inflate. Hence, we can choose accordingly the right 
linearized solution^. Taking the above arguments into account, we now turn to the 
perturbative considerations. 

Expanding the metric around the 5D Minkowski vacuum 

(5) _ | 7 

9mn — VMN -r tlMN , 

and using the harmonic gauge in the bulk, we obtain that the only nonzero compo- 
nents, h,j, v and /155, satisfy the relation 

h = K = h\ , 

where the brane is located at 25 = 0. The linearized Einstein's equations take the 
form p] 

^ D 5 + ^ 8(x 5 ) D 4 ) V = (t^ - iiv) + ^p M cU , (2.3) 

while the trace of this equation is 

M 3 1 

-±n 5 h = -6{x s )T. (2.4) 

For a domain wall T = — 3a 8 (z). By the symmetries of the problem, the static 
solution should depend only on |x 5 | and \z\, and satisfy d 5 h = for x 5 = (and 
z 7^ 0). One can show that up to a constant, the only such a solution with the 
appropriate singular behavior at z = x 5 = is 

h = -l^lo g [(z 2 + 4)^} , (2.5) 

where \l is a mass scale, which when the transverse size d is incorporated turns out 
to be of order 1/d. Once we found h, we can plug it into the r.h.s. of (I2.3p . to 



6 Note that we would obtain a wrong linearized solution if we did not know the exact solution 
and we were to follow the intuition gained on domain walls in 4D GR. 
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find the full h^ v . Now, since on the brane h depends only on z, one immediately 
sees that the components parallel to the DW have no source term. Requiring that 
they are regular at infinity, one obtains that they should vanish. Thus, h = h z z , and 
indeed the (zz) equation is trivially satisfied. Hence, to leading order in a /Ml, the 
metric is 

ds 2 = (1 + h)(dx 2 5 + dz 2 ) - dt 2 + dx 2 + dy 2 , (2.6) 
with h given in (12.51) . The metric (I2.6P can be brought to a locally flat form 

ds 2 = dp 2 + p 2 d6 2 - dt 2 + dx 2 + dy 2 , (2.7) 

with the range of the angular coordinate modified (to leading order) as < 9 < 2tt—5 
with 

6 = Wy 

This is the well known conical space induced by a local codimension-2 object, with 
deficit angle given by 5. Interestingly, this feature holds true also in the exact 
solutions for sub-critical domain walls (those with o < Ml) on the conventional 
branch, as we shall see in Sec 13.11 Let us remark now that (I2.6P implies that 
the induced metric on the brane is flat; setting x§ = and using the coordinate 
dz* = a/1 + hdz, the induced metric is dz 2 — dt 2 + dx 2 + dy 2 . Hence, a (sub-critical) 
domain wall in DGP has no gravitational effect on the brane. As will become more 
clear below, the wall switches on the extrinsic curvature term in such a way that 
the latter behaves as a negative tension wall, screening it completely. 

3 Exact solutions 

The DGP action pQ of the bulk plus brane plus the Domain Wall (DW) in the thin 
wall approximation takes the form: 




where R$ and R4 are the Ricci scalars of the metric in the bulk g^J and of the induced 
metric on the brane g^, and 7^ is the induced metric on the DW. We have also 
introduced the brane tension r, even though our main interest is in the case r = 0. 
In the derivation of the exact solutions below, we follow [15]. Given its potential 
relevance for the present acceleration of the universe [HI [17], we shall discuss both 
the conventional and the self-accelerated branches of the theory [HI [T9] . 

The equations of motion arising from (13.11) can be split into the equations for the 
bulk and those for the brane. We assume that there is no cosmological constant in 
the bulk, and that it has at least the same symmetries as the (maximally symmetric) 
domain wall. We are mostly interested in the cases when (1) the DW is inflating, 
in which case its geometry is a 3 dimensional de Sitter space; (2) the DW is flat. 
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Then, the analog of Birkhoff's theorem in 5D ensures that the bulk metric locally 
has the form 

ds 2 = f{R)dZ 2 + ^ + R'dsl , (3.2) 

with f(R) = u — C u /R 2 and C u an integration constant. Here, ds 2 is the line 
element of a three dimensional Minkowski space (u = 0) or de Sitter space with 
unit curvature radius (u = 1). In this article, we shall consider only those boundary 
conditions in the bulk that fix it to be flat, that is u — 1, C — 0. The general case 
will be discussed elsewhere [20]. Still, it is convenient to introduce two different 
slicings, 

ds 2 = dZ 2 + dR 2 + R 2K ds 2 K , (3.3) 

corresponding to the full Minkowski space (k = 0) and to a kind of Rindler space 
(k, = 1) . As we will see shortly, the 3-dimensional geometries of the domain wall 
are given by ds 2 . Hence, the Rindler (Cartesian) coordinates are well suited for an 
inflating (flat) wall. 

The equations for the brane are given by the Israel junction conditions, 

2M*K, u = Tf u -±T cS g, u , (3.4) 

where we imposed Z2 symmetry across the brane, K^ v is the extrinsic curvature, 
T efr = T eS^ and 

Tf u = -m 2 P G$ - rg, u - a5(0l»» • (3-5) 

Here, £ is the coordinate along the brane orthogonal to the DW, and r is the brane 
tension. 

The idea is now the following. Once the form of the bulk is known, the problem 
reduces to finding the embedding of the brane in this space, via Eq. (13.41) . The 
location of the brane can be parameterized by two functions (R(^),Z(^)). Then, 
the induced metric on the brane is 

dsl = de + R 2K {i)ds 2 K1 (3.6) 

where we have chosen the 'gauge' 

Z 12 + R' 2 = 1 , (3.7) 

and ' denotes a ^-derivative. This can be thought of as fixing the arbitrariness intro- 
duced in the parameterization (R(^),Z(^)). We then compute Eqs. (13.41) explicitly 
in terms of this parameterization, which allows us to solve for -R(£) (and Z(£), 
through (13.71) ). Generically, the brane will 'cut' the Z — R plane in two regions, and 
the full 5D space is then obtained by taking two identical copies of the appropriate 
side pasted together along the brane. 

The components of (13.41) along the 3D de Sitter/Minkowski directions lead to 



VI - R' 2 l-R' 2 _ 



eK ^-R-= -"^^ + «' (3 ' 8) 
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where 

r.= if ■ (3-9) 

The parameter e = ±1 is the sign of the extrinsic curvature. For r > 0, e is positive 
(negative) for the Conventional (Self- Accelerated) branch. 
Equation (13.81) leads to 

1 _ D/2 

~^-^H\ (3.10) 

where H solves 

e k — — = — k -\- — — =- . (3.11) 

r c 3mp 

The case k = 0, can only be met for r = 0, which makes perfect sense, since in this 
case the metric (13.61) is flat. For k = 1, the solutions symmetric across the DW are 
of the form 

R(0 = H- 1 sm[H(£ -m (3.12) 
where we have introduced the integration constant £ so that the DW is placed at 

e=o. 

The (££) component of fj3.4j) . involving the localized source term ~ cr5(£) deter- 
mines the constant £ - This equation is not independent of Eq. (13. 8p . which is valid 
away from the source localized at £ = 0. Close to the source we obtain the following 
divergent terms 

TD/I Till O 

- c ^= = ^ + —m , (3-13) 

where 

r dw = ^ , (3.14) 
a 

is the would-be domain wall horizon. Integrating both sides of H3. 131) from £ = — e 
to £ = e, and letting e — > 0, one obtains the following matching condition for the 
jump of i?'Jll 

e /\R f 2 
Aarcsini2n = K H , (3.15) 

2r c R r dw 

where AX = lim e ^ [A(£ + e) — X(£ — e)]. To arrive at (I3.15p . we have used that i? 
is continuous at £ = 0, and it can be treated as a constant for £ close enough to 0. 

Now, the gauge condition (I3.7p can interpreted as the condition that the vector 
(R',Z') has unit norm. We can thus identify \R' \ = sin/5. From Fig. [TJ /3 is the 
angle between the tangent to the brane and the Z axis at the DW location. Upon 
cutting and pasting the appropriate sections of the R — Z plane, we realize that (3 
is related to the deficit angle (see Fig CD) through 

5 = 4e(3 . (3.16) 



7 Note that in Ref. [15] the integration of the l.h.s. of (|3.15[) is incorrectly done. 
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P = tt/2 - a 



BULK (SAB) 




BULK (SAB) 



Figure 1: The brane trajectory in the R — Z plane, for a generic DW in an inflating 
brane (H ^ 0). The two figures represent a constant time snapshot, and one of 
the directions along the DW (represented by the thick dot) is obtained in the upper 
diagram by rotating it around the Z (vertical) axis. The space is cut by the two 
constant Z sections where the brane is at an angle a = H^q. The space in between is 
then removed, and the two slices joined together. Then, another identical copy of the 
resulting 'pancake' is glued along the brane, as indicated in the lower diagram. For 
the conventional branch (CB), the bulk is the interior, while for the self-accelerated 
branch (SAB), it is the exterior. For the CB, this space has a deficit angle 5 = 4/3, 
with (3 = 7r/2 — a. For the SAB, we have an excess angle, that is 6 = —4/?. 



Note that the sign of 5 tells us whether it is a deficit or excess angle. Equation 
( 13.15P then takes a simple formal 

Rq 4:r c r dw 

The interpretation of this is clear. It tells us how much the DW is behaving as a 
codimension-2 object (opening-up a deficit angle), versus a codimension-1 object 
embedded on the brane worldvolume (generating a jump in the derivative of the 
induced metric). 

Assuming Z 2 symmetry across the wall Eq (13.171) reads, in terms of (3, 

2k Hr c tan 3 + e (3 = . (3.18) 

4M| v ' 

This equation cannot be solved analytically (except for k = 0), but it is straight- 
forward to solve it numerically for given values of the parameters. We shall now do 



8 This correctly reproduces the 4-dimensional and 5-dimensional limits. For r c — ► oo, we obtain 
—AR' /Ro = 2 /rd w > which we identify as the Israel junction condition for a DW in GR. For r c = 0, 
we get S = liirL-.— »o 8r c /r,j, w — <r/M^ , which is the correct relation between the tension and the 
deficit angle. 
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this in the next subsections. 

For clarity, we shall now summarize the main features of these solutions. It is 
already apparent from (13. 17j) that both the 5D effect (amount of deficit angle) and 
the 4D effect (the jump within the brane induced metric) produced by the wall will 
in general differ from the usual pure 4D or 5D ones. In order to quantify these 
deviations we can always introduce the effective 4D and 5D tensions in terms of 
whichever geometrical effect is produced, 



and 



(ad) 2 A-^o 
cr ff ; = —2k m P 

~o 



(5D) 



(3.19) 



a^-'=5Ml. (3.20) 

Recalling that the deficit angle in (13. 1 T[) arises from the extrinsic curvature and that 
the jump in R' comes from the intrinsic curvature, we readily identify that <J^ D ' 5D ^ 
contain a contribution from K^ v and from respectively, aside from the DW 
tension itself. Also, in terms of cr^ ' 5 , the wall junction condition (13.171) can be 
expressed as a 'sum rule', 

+ = ° ■ (3-2i) 



(AD) 
eff 



Note also that the notion of the 5D tension (13.201) is not valid for super-critical walls 
(a > 2ttM%), simply because the deficit angle cannot exceed 2n without adding 
extra sources. The appropriate notion in this case will require some previous famil- 
iarization with super-massive codimension-2 objects, and will be discussed at length 
in Section [3731 



branch 


DW tension 


(AD) i 


(5D) 1 


CB 


a < 2-nMl 





1 




a > 2ixMl 


, 2ttM^ 

1 - 

a 


d 

2^c 


SAB 


a < 2txMI 


2(l- (a/M, 3 ) 2 /48 + ...) 


-(l-{a/M!) 2 /24 + ...) 




a > 2-nMl 


2itMt 
1 + - + ... 

a 


d 

~2^ c 



Table 1: Summary of the screening properties of our solutions. See the discussion 
in the text for the definitions of (J^^ and cr^ . Notice that for sub-critical walls 
the sum rule (I3.2ip is satisfied. 



Now, let us briefly go through the results displayed in Table [TJ As mentioned 
in the introduction, in the Conventional Branch, sub-critical walls suffer a complete 
screening of their 4D tension. The wall effectively switches on the extrinsic curvature 
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in such a way that it completely compensates for the wall tension. Let us emphasize 
that the 4D curvature of the brane is zero in this case even on the DW location. 
Precisely for this reason, the 5D tension is not screened at all. 

Still for the Conventional Branch, the supercritical walls show a very different 
behaviour. First of all, the wall starts to inflate. So, the screening of the 4D tension 
is not total. At the same time, the 5D tension is screened by a huge factor. What 
happens in this case is that the space transverse to the wall is compactified to a 
size of the order of the thickness d of the DW core (we are assuming that d is 
finite but small, d -C r c ). Hence, the model changes significantly, e.g., there exists 
a graviton zero mode. Nevertheless, as we shall show in Section 13.31 'ordinary' 
super-massive codimension-2 objects (i.e., without the DGP brane) also compactify 
one direction of the transverse space, and start to inflate. Hence, the 5D effect to 
which we should compare is the worldvolume inflation rate, rather than the deficit 
angle (which is saturated to 2tt) . As we shall see, the walls in a DGP brane inflate 
at a rate suppressed by d/2r c relative to the ones without the brane. As argued 
before, the analog of r* scale for for all kinds of walls coincides with the thickness 
d. Hence, quite surprisingly, the super-critical walls display the same screening of 
the 5D tension/mass found in [7] for the Schwarzchild-like objects. 

In the SA branch, we will see that the deficit angle is negative, meaning that 
the 5D tension is 'over-screened' (the sign of the net tension a^ D > is reversed). On 
the other hand, from the 4D point of view, there is anti-screening: the extrinsic 
curvature takes the form of a positive tension wall. For small tensions, mimics 
a wall with almost exactly the same tension, so the 4D effect is enhanced by a 
factor 2. In the super-critical case, we shall find a suppression of the Hubble rate 
on the wall ~ d/2r c , while the deficit angle is close to — 2n instead of 2tc (which is 
what the sign in the lower-right entry of Table [1] stands for). 

3.1 Conventional Branch 

This branch is defined by requiring that the brane extrinsic curvature be positive 
(e = 1) for non-negative brane tension r. As seen from ( I3.16p . this implies that a 
positive tension wall generates as usual a deficit angle in the bulk. From (13.171) . 
we also see that this can be interpreted from the point of view of the brane as a 
screening due to the extrinsic curvature term, which behaves like a negative tension 
domain wall. As we will now see, the amount of screening depends on the brane 
tension, and in the extreme case r = 0, there is complete screening. 

Flat domain wall (k = 0) 

In this case, Eq. (13.81) is automatically satisfied if we set r = 0. On the other 
hand, Eq. (I3.13P is of the form R" oc Hence, taking into account (13.71) . the 

trajectory is of the form 

fl(0 = |£l sin/3 Z(£)=£ cos/3 (3.22) 
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The junction condition on the DW (13.171) fixes /3, and leads to an expression for the 
deficit angle 

6 = Wy 

as for an ordinary codimension-2 object. This space is illustrated in Fig. El It makes 
sense as long as the deficit angle is less than 27r, that is for 

a < 2nMl . 

We recognize that the linearized solution of Sec. [2] was tailored to reproduce this 
exact metric. 




Figure 2: Conical space corresponding to the flat DW case in a flat brane. The 
deficit angle is 5 = 4/3, which is always less than 2tt. This figure refers to the 
Conventional Branch, and with zero brane tension. 

Now that we have a geometrical method to find the solutions, we could easily 
construct multi-center solutions with as many DWs as we'd like. The induced metric 
is always trivial (up to the topology), and the embedding of the brane in the bulk 
has a kink-shape, as in Fig [2j at each DW. In particular, if the wall tensions add 
up to 27r, then the bulk is compactified with compactification radius related to the 
separation between the walls. 

Let us now describe in more detail the critical case, o = 2nM^. A single 
critical wall leads to a seemingly singular space because the angular component of 
the metric (12. 7\\ vanishes. This is of course an artifact of the thin wall approximation. 
If we introduce the thickness of the wall d, then the metric asymptotes that of a 
cylinder of radius d. Hence, the bulk (the transverse space) is compactified, and in 
the 'thin wall' limit d — > 0, one transverse direction disappears. We can illustrate 
these statements giving the critical wall some structure, for example by splitting it 
into two walls with tensions 37r/2 and tt/2 (in units of Mf), separated by a distance 
d. Since each of them is subcritical, the bulk is a wedge of 2n minus the deficit angle 
at each wall. In this case, instead of Fig the bulk looks like (two copies of) a 
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semi infinite strip of width d. Thus, the compactification radius coincides with the 
separation between wallgjl 

This solution contains one modulus, d. At first, this is puzzling because it seems 
that it costs no energy to move the walls apart, but this would result in changing 
the compactification scale. The reason why there is no paradox is that to change 
globally the position of the walls actually costs an infinite amount of energy. The 
fact that d is a flat direction means that local perturbations in the wall separation 
propagate as massless fields. 

There appears to be another paradox, though: the bulk is compactified, so one 
expects to have four dimensional gravity at distances larger than d. However, the 
DW is flat, so it does not gravitate as in 4D. This situation occurs also in ordinary 
Kaluza-Klein theories, without the brane. A critical codimension-2 object in 5D 
compactifies one transverse in the same way, producing an effectively 4D space. 
But it does not inflate, as it would in 4D. This an example of off-loading: all the 
gravitational effect of the DW tension is 'exhausted' in compactifying the transverse 
space. We will see below that if the tension is supercritical, then it starts behaving 
as in 4D, and in the limit a ^> M% we will recover the usual metric of a domain 
wall on the brane. But in the critical case, the four dimensional behavior is not 
reproduced. 

Finally, we shall mention that the space generated by a critical (distribution 
of) walls is unstable under nucleation of bubbles of nothing, very much like with 
the Kaluza-Klein vacuum [21]. The reason is that, the space is asymptotically 
cylindrical, which has the same topology as the bubble of nothing instanton. The 
estimated nucleation rate is of order exp(— m 2 P R 2 ) [21], where R is the (asymptotic) 
radius of the extra dimension. In our case, this is related to the wall thickness. For 
a single critical wall, assuming that its thickness is of the same order as the scale of 
the tension M*, this gives a highly suppressed rate. For an extended distribution of 
walls R is larger, so the rate is even more suppressed. 

Inflating domain wall (k — 1) 

Solving for f3 in Eq. (I3.18P for r ^ 0, we can find the deficit angle 5 and the DW 
radius Rq = H^ 1 cos (3. The results are shown in Fig. [3] for Hr c =1/5 and 1/100. 
We can see how as H — > 0, the subcritical walls (with a < 2ttM^) are much flatter 
than in usual GR. 

For H = 0, Eq. (13.101) implies that R' 2 = 1. Thus, the brane trajectory lies 
along Z = const, lines in the Z — R plane, which means that the deficit angle is 2ir 
(for the CB). Strictly speaking, this limit is singular, because the bulk collapses to 
zero volume. However, if we introduced the DW thickness d, we would realize that 
the bulk is a slice of size roughly given by d. Here, we shall ignore these details and 
assume that the geometry is appropriately regularized in such a way that it still 

9 One can see that for tensions giving individually deficit angles Si < tt and 2ir — Si the com- 
pactification radius is dsin<5i. 
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Figure 3: The radius of the DW Rq for the conventional branch (CB), in the Z2 
symmetric case, for Hr c = 0.2, 0.01 and 0, and the deficit angle (for Hr c = 0.01) 
as a function of the tension a. We see how the DW tension is screened by the K^ v 
term (contributing effectively negative tension), so that the radius of the DW is 
larger than in 4D GR (dashed line). There is a qualitative change of behavior for 
tensions below or above the critical value 27rM^ (vertical lines). In the limit H = 0, 
tensions smaller than or equal to 2irM^ induce a flat DW. Only for o > 2irM^, the 
DW worldsheet is inflating, and Rq approaches the 4D value (raw) for o M% . 



makes sense to speak of a deficit angle 2tt. Then, the matching condition (I3.17P 
takes the following form, 

Ro 4r c r dw 

The solutions of this equation are represented in the lower left plot of FigO Solutions 
with a finite Hubble radius Rq start to exist for supercritical walls o > 2irM%. 
Subcritical walls have an infinite Rq, that is they are flat, in agreement with the 
arguments of the previous sections. 



3.2 Self-accelerated branch 

In this case, the brane extrinsic curvature is negative (e = —1). This implies that a 
positive tension wall generates a negative deficit angle in the bulk (see FigfU or Eq. 
( I3.16P ). From the point of view of the brane, this is seen as an anti-screening due 
to the extrinsic curvature term. Using e.g. ( 13. 18[) for tan/? ~ (3 <C 1, one obtains 
tan/3 oc R' /Rq twice as large as in GR. This means that in this case K^ u behaves 
like a positive tension domain wall with effective tension a. 
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From the point of view of the bulk, however, there is 'over-screening'. Indeed, 
the deficit angle is negative (see FigH]), which means that the effective codimension- 
2-tension is negative. This can also be visualized from FigdJ for the SA branch, the 
bulk is the exterior of the brane, which gives an excess angle. 

We can also see this from Eq (13.181) for small DW tension, o <C M*. In this limit 
also /3 <C 1, so the deficit angle is approximately (here we take k — 1) 

1 a (3.24) 



1 + 2eHr c Ml ' 

Hence, in the CB there is screening, while on in the SAB there is 'over-screening'. 
In fact for the SAB with zero tension, the deficit angle is exactly opposite as the 
naively expected. 

Two comments are now in order. Given that from the point of view of the 
bulk the wall in the SAB behaves as with negative tension, it appears that the 
stability of this solution is not granted. However, a perturbative analysis would not 
be conclusive, since the brane is in a strong coupling regime for this solution [221 EES] ■ 
On the other hand, it is clear that a perturbative treatment as in Section [2] will not 
capture the over-screening. This is precisely a consequence of the non-perturbative 
nature of the solution. 




-61 -10 -5 5 10 4 M 3 -10 -5 5 10 4 jy[ 3 



Figure 4: The deficit angle and the radius Rq of the DW for the Self-Accelerated 
Branch. The radius i?o approaches its 4D value (dashed line) for large a. The 
fact that Ro is smaller than the GR result is a consequence of the anti-screening 
produced by the extrinsic curvature term. In the third plot, we see the anti-screened 
the 4D tension, which amounts to a factor 2 for small a. For large tension, a^ 4D ^ 
approaches the GR result. 



3.3 Super-critical walls and screening of the 5D tension 

In this section we discuss the spacetime produced by a (local) supermassive co- 
dimension-2 brane in ordinary gravity, in order to compare it to the supercritical 
walls found above. Solutions of a super-massive Abelian-Higgs vortex have been 
previously found numerically [231 121] • The metric was shown to be regular and 
time-dependent, essentially because topological- inflation sets in. 
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In our case, the codimension-2 object is assumed to have the form of (thick) a 
Domain Wall confined on a 'fake' brane (with negligible tension and thickness, and 
no induced gravity term). As such, it does not have cylindrical symmetry. Rather, 
it is like a strip of a codimension-1 object of a finite extent, which, seen from 
afar appears as codimension-2. Thus, the problem under consideration is slightly 
different from that of [231 1211 • The advantage is that we can find the solution 
analytically, and the detailed structure of the core is not going to be important. 
Let us also emphasize that the solutions presented in this section are not the most 
general compatible with the assumed symmetries, since we will also require that the 
space outside the object is locally flat (this is a similar to the restriction to C u = 0, 
u — 1 in (13.21) ). We shall show that also these ordinary codimension-2 objects with 
tension a > 2ttM^, the worldsheet starts to inflate, and we shall compute at which 
rate. 

We can proceed in the same way as in Section [3] with vanishing brane tension 
t = 0. We shall concentrate on the five dimensional case, and take a generic form 
of the energy momentum tensor compatible with the symmetries of the wall, 

Tjl = diag(-p, -p, -p, P) , 

and we assume that p and P (if nonzero) depend on £ only and are peaked around 
£ = 0. We shall not assume any particular microscopic model for the wall. Let us 
just mention that thick gravitating domain walls arising from scalar field models do 
develop a non- vanishing pressure in the core [25] (see also [26]). 

The embedding equations (I3.8P and (13 . 13[) without the induced gravity term and 
with t = become 

2. Ml = , (3.25) 

and 

which in fact is equivalent to (I3.25P once the conservation equation 

P' + 3k ^ (P + p) = 
R 

is used. 

Equation (13.251) immediately tells us that P < 0, and that it should vanish when 
the worldvolume does not inflate [k = 0). Instead, if it inflates, then P is nonzero 
essentially at the location of the object, where R' 2 ^ 1. 

As seen from (13.221) . a subcritical wall produces a discontinuity |Ai2'| < 2 while a 
critical wall saturates this bound. For supercritical walls, then there are two options: 
one is to have a larger |AJ?'|, but this would require another source at finite distance. 
The other is that the P switches on. Let us now see that this is compatible with 
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keeping |AP'| = 2 for the supercritical walls. The integral of Eq. (I3.26P around the 
core (that is from B! = 1 to R' = — 1) gives 



2vMl = J d£(p + 2P/3) 



If we keep calling a = J core p <i£ the tension, and defining the width of the wall d by 
f Pdti = Pn d, we obtain 



^P d = a -2ttM! 



3 

where P = P(0). Now, Eq. (13.251) tells us that at the center of the core (R' = 0) 
the Hubble rate 1/Po is related to P as 

— Pn 



Po 

we arrive at 

i- = —I— ( a - 2tt Mf) . (3.27) 

Given that P' 2 = 1 outside the walls the transverse space is compactified to a length 
equal to the thickness d. Hence, the theory becomes effectively 4-dimensional, with 
an effective 4D Planck mass 

Ml = Mld . 
With this in mind, we can rewrite (I3.27P as 

R~o = AM! ~ 2d ' (3 ' 28) 

Hence, we see that we obtained the expected results: Eq (I3.27P tells us that for 
tension larger than 2irM^, the wall starts to inflate; The Hubble rate approaches 
that of a codimension-1 object in GR (in terms of the effective Planck mass M4). 
The interpretation of (13. 28ft in 4D terms is also interesting: the 4D tension of the 
wall is screened by the higher dimensional effects by an effective mirror wall with 
negative tension and an equivalent horizon of order d. 

It is also worth describing the structure of the space-time generated by this ob- 
ject. The space is (the interior of) a flat 'pancake' of a thickness essentially given 
by the thickness of the wall, d. Hence, the picture is similar to Fig. (TJ but with 
the 'sides' completely flat. The finite expansion rate on the wall translates into the 
fact that the object is placed at a finite radial distance P (recall that in FigHJ the 
horizontal coordinate is a radial Rindler coordinate). Accordingly, the presence of a 
horizon at a finite distance from the object is seen in the picture because the coor- 
dinate orthogonal to the object 'terminates' at the origin P = (which is actually 
the horizon). 
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Now that we know how supermassive codimension-2 objects behave, we can 
compare to what we found for the DGP model. Note that (I3.28P is very similar to 
(I3.23P with the thickness d instead of 2r c and M 4 instead of mp (the coefficient in 
the induced gravity term). One immediately realizes that the Hubble rates induced 
on the objects in DGP is extremely suppressed as compared to ordinary 5D gravity, 

~ l \nT = — ■ (3-29) 

Hence, the 5D tension is suppressed by a factor d/r c , which as argued before is 
the same as r*/r c , since the scale where the nonlinearities build up is r* = d for 
domain walls. Note that this is exactly the same suppression that occurred in the 
Schwarzchild solution found in [7], which also applies to the supercritical case where 
the radius of the source is less than its own r*. Note also that for the domain walls, 
the 4D tension is screened, even though not completely. Indeed, without screening 
of the 4D tension one would expect 1/Rq = a /Amp and instead one has (see (I3.23P ) 

1 a 7i 



Ro Amp Ar c 

At this point, it might seem that the screening of the 5D tension (13.291) for super- 
critical walls is suspiciously discontinuous: in the 5D limit, r c — > 0, (13.291) should 
go to 1. Let us now show that this is precisely what happens, once we include the 
thickness of the wall from the beginning when the wall is on a DGP brane. We only 
need to repeat the above steps keeping the induced gravity terms. The junction 
equations are 

„ r , Vl - R' 2 o 1 - R' 2 P , 

and 

On one hand, we can read from (I3.30p that 



+ ^ • (3-32) 



_Po = ( 2eMl | m 2 P 
3 \ Rq i?o 

Integrating this equation along the core of the wall, we obtain, for the Conventional 
Branch, 

2-nMl = ldP + a - Am 2 P ^- . (3.33) 

From this and (I3.32p . we arrive at 

m p + d Ml m 2 P d , , 

4-^ *- + 2-|p = a - 2tt Ml . (3.34) 
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Assuming that d <C r c , one sees that the for Rq 3> d the 1/Rq term dominates. 
Hence, we obtain 

1 1 a-2nM* 

W ~4m 2 P + dMf (3 ' 35) 
which gives the suppression factor 

d d 



d + 2r c 2r c 

Note that (13.351) being consistent with our assumption that Rq 3> d leads to 



(3.36) 



8 1 , 
d < - - r c . (3.37) 

This, though, does not seem to be much of a constraint for microscopic walls unless 
the tension is extremely large. 
For the SAB, fl3~33l reads 

- 4 arccos(i2g/r^) Ml = \d P + a + 4mp — (3.38) 

3 R 

so the analog of ( 13.341) 

4mp 1 ~5 2rC + 2^ = a + 4M, 3 arccos(^/r c 2 ) . (3.39) 

Rq Rq 

Again, assuming that d Cr c one obtains for d i?o "C r c , 

1 a + litMl 



Rq 4(1 - d/2r c )mp ' 
leading to the suppression factor 

a + 2nM? d 



(3.40) 



a - 2-nM'l 2r c - d 



which, for a ^> Ml coincides with the one for the CB, (13.361) . Also, the consistency 
of (13.401) with the approximation d Rq r c leads to a mild condition similar to 

(E3ZD. 



4 Domain walls in IR modified gravity 

Our findings of the previous section may well hold in more general models of modified 
gravity. Detailed arguments in favor of this will be discussed below. However, before 
we plunge into those considerations we'd like to emphasize the following: Nonlinear 
effects played an important role in the previous section where we obtained the results 
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in DGP. Unfortunately, it is not clear whether the general models which we will 
discuss below have any consistent nonlinear counterparts. For one, massive gravity 
is known not to have a consistent nonlinear completion [27J EH1 EH] , as any of its 
versions gives rise to certain ghost-like instabilities. All of our discussions concerning 
general models below ignore these difficulties, and we have no way of assessing at 
present the significance of this for our discussions. 

The most general ghost-free linearized theory in which gravity is mediated by a 
symmetric tensor field h^, is a generalization of the Fierz-Pauli [30J massive theory, 
and satisfies the following equation [13], 



£^Kp - m 2 (D) (V - n^h) = -16nG N , (4.1) 



where, 

e$h af) = DV - n V»»h - d a d,h au - d a d v h ait + ri»ud a d?h a p + d^h , (4.2) 

h = and Gn stands for the Newton's constant. 

For m 2 (D) = this reduces to the linearized Einstein's equations, describing the 
propagation of two polarizations of a massless spin-2 graviton. Here, we shall be 
interested mostly in the case m 2 (D) = const, that corresponds to massive gravity; 
we will compare the latter to m 2 (D) = ^/0/r C7 which corresponds to the DGP 
model. Our conclusions also apply to the more general cases where 



m 



2 (D) =r; 2{1 - a) D a , (4.3) 



with < a < 1 (the lower bound arising from unitarity [13J), assuming that these 
models have sensible nonlinear completions. 

Gravitons satisfying (14. II) propagate 5 polarizations for nonzero m 2 (D), which is 
at the root of the vDVZ discontinuity [2j. The naive linearized metric produced by 
a localized source T Mi , is given by 

V = - 16 ^ D _1 2(D) - I (v - -^o) 9 ^) T ] > ( 4 - 4 ) 

where T = Tff. The factor 1/3 reflects the existence of the extra helicity-0 polariza- 
tion x, residing in h^ u . The latter could be separated as follows: 

V = V + \nnv X + I J^y X , (4-5) 

where contains the usual 2 polarizations, as in the massless case, and the helicity- 
1 states are ignored as they do not couple to conserved sources in the linearized 
theory. Convoluting Eq. f)4.4p with a stress-tensor of a test source T' , we obtain 
the following one-graviton exchange amplitude 
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Let us apply this to an infinite domain wall with the stress-tensor given in (11.11) . 
Taking a non-relativistic probe T' = M 5°<5° S^(r), we obtain that the amplitude 
equals to zero! Actually, it is also easy to see that it vanishes for any conserved 
source. First, note that for the DW, the tensorial structure T^ v — (l/3)T?7 Atl/ oc 
Hence, the amplitude is 

J d 4 x T' z J{z) = - j d'x d z T' zz j Z dz'f(z') = J d*x d A T' Az J* dz'f(z') = , 

where f(z) is the result of applying (□ — m 2 (D))~ 1 on 5(z) (certainly a function of z 
only). In the derivation above, we integrated by parts, used the local conservation of 
T' , and integrated by parts again (the label A stands for all coordinates except z). 
Hence, we conclude that the amplitude vanishes, which is due to the 'pure gauge' 
form of the metric (I4.4|) for the DWs. 

Thus, in the naive linearized approximation, domain walls in massive gravity 
(and in general theories parameterized by (14.31) ) do not gravitate. This conclusion, 
however, is not warranted until the nonlinear effects are taken into account. This 
is because of two interrelated reasons already discussed in Sections 1 and 2: (1) 
There is an ambiguity in linearized solutions; one cannot decide whether to choose 
a time-dependent or a static solution, without knowing what is going on with the 
source itself. (2) In modified gravity, non-linearities introduce a scale r* at which 
the perturbative expansion breaks down, and, therefore at r <^ r* the above results 
cannot be applied. In DGP these issues were addressed by finding exact nonlinear 
solutions in Section [31 We do not have the luxury of a stable nonlinear theory in 
the case of massive gravity (and the general a-theories (14.31) ). However, this may 
only be a temporary technical problem, and if so, it would make sense to try and 
estimate what the analog of the r* scale for a domain wall in massive gravity is. We 
will do this below. 

First we recall that the r* scale and strong coupling originates from the l/m 2 (D) 
term in (14. 4p . which at the non-linear level gives rise to interactions that become 
singular for m 2 (D) — > |3j. In a schematic way, the leading singular vertex has the 
form [13] 

rt {1 - a) ^x(d X )\ (4.6) 

where we should keep in mind that the four derivatives do not necessarily come in 
the form of D 2 . For a = 0, 1/2, this was found in [3] (see also [31] for a = 0). Due to 
this singularity, the localized gravitating sources end up having a new 'horizon' at 
r = r*, where the perturbative expansion in Gn breaks down. Consider now a non- 
relativistic spherically symmetric source T^ u = M 5®5® 5^ 3 '(r), with a gravitational 
radius r g = IG^M. The scale r* is then found as follows: Assume that far enough 
from the source, the linearized solution 

_ i _ gggjj 
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is valid. From this and (14.51) . one finds that for r Cr c 

X=4^- (4-8) 
2 r 

At the scale r* the contributions from the nonlinear diagrams (14. 6 j) catch-up with 
the leading linear one. This happens when 

r c \4(i-«)r g 
— I — ~ 1 . 
r / r 

This determines the Vainshtein scale [32] 

r v ~ (r 9 m- 4 ) 1/5 , (4.9) 
for massive gravity, and the r* scale 

^( Vc 2 ) 1/3 , (4.10) 

for DGP. At r ^ r*, the metric can be found as a series in r/r c [1]. 

Let us now apply similar considerations to subcritical domain walls. Assuming 
that the linearized solution 

5 Z 5 Z + 9ttdu 



Ku = 3 2 7n? Wk<tG n 8{z) , (4.11) 

□ — m z [U) 



is valid far-away from the source, for z <C r c Eq. (14. lip leads to 

X = — 87T GVcr |^|. 

With such a dependence on z the nonlinearities vanish due to the structure of the 
vertices. This is obvious for the leading singular vertex (14.60 . and is also true for 
all the other interactions. The key point is that both for m 2 = const and for 
m 2 (D) = \/TJ/r c , all non-linear vertices include at least one x ( or fy*") with at 
least two derivatives. As a result, for solutions for which h^ v is linear in z, the 
strongly coupled vertices vanish outside of the source. The only place where they 
can make contributions is in the core of the source itself. Thus, the scale r*, for a 
subcritical domain wall, is determined by the transverse size of the wall itself. This 
is reminiscent of a metric of a subcritical Schwarschild source [7], except that the 
subcriticality condition in the two cases are quantitatively different. 

Finally, let us discuss a similar issue for a wall with a finite longitudinal extent. 
Specifically, consider a disk of radius L located on the plane z = 0, and let us measure 
the gravitational interaction of this 'wall' along the z axis. For r, C z <C I, the 
solution should be well approximated by (14. lip , and x \ z \ as before. In this case, 
we could apply the same reasoning for strong coupling vertices as we did for an 
infinite wall, and could conclude that they vanish outside the wall as long as r* is 
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smaller than L. Let us get a naive estimate for the value of r* for a finite wall. The 
total 'mass' of the source is (here we ignore the fact that the domain wall also has 
a negative pressure) 

M ~ naL 2 ~ it5M^L 2 , 

where we introduced 5 = a /Ml, and M* is the fundamental Planck scale in the bulk 
(let us concentrate on the DGP model for the moment). Then the naive estimate 
r naive _ ^G^Mr*) 1 / 3 would suggest 

r naive _ 1/3 _ 

If r" arae is smaller than L, then non-linearities never contribute outside the wall. 
This gives the following constraint on the size of the wal0 

L>5r c . (4.12) 

Hence, in order to avoid the strong coupling of the longitudinal mode x> the wall 
radius has to be quite large; for 5 ~ 1 it should be larger than r c . This is consistent 
with the expectation that the behavior of a DW here is different from that in GR. 
The above follows from the fact that the infinite domain wall is already probing the 
scale of modification of gravity r c . 
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